Suppression of Anderson localization of light and Brewster anomalies in disordered 
superlattices containing a dispersive metamaterial 
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Light propagation tlirougli ID disordered structures composed of alternating layers, with random 
thicknesses, of air and a dispersive metamaterial is theoretically investigated. Both normal and 
oblique incidences are considered. By means of numerical simulations and an analytical theory, we 
have established that Anderson locaUzation of light may be suppressed: (i) in the long wavelength 
limit, for a finite angle of incidence which depends on the parameters of the dispersive metamaterial; 
(ii) for isolated frequencies and for specific angles of incidence, corresponding to Brewster anomalies 
in both positive- and negative-refraction regimes of the dispersive metamaterial. These results 
suggest that Anderson localization of light could be explored to control and tune light propagation 
in disordered metamaterials. 

PACS numbers: 42.25.Dd, 78.67.Pt, 72.15.Rn 



CZ3 



I 

o 



> 

o 
o 

o 
o 



X 



In the last decade or so, advances in metamaterials 
have allowed the development of unusual optical proper- 
ties, with no counterpart in natural media, opening up 
new frontiers in photonics. For instance, metamaterials 
exhibit negative refraction [ij , resolve images beyond the 
diffraction limit Q , exhibit optical magnetism 0, 01 j act 
as an electromagnetic cloak 0, Q and yield slow light 
propagation 0. 

Notwithstanding the wide range of new physical phe- 
nomena unveiled sofar, there are fundamental issues that 
are still not fully understood. Anderson localization (AL) 
of light in disordered metamaterials is certainly one of 
them. The concept of AL was originally conceived in the 
context of condensed matter physics as the vanishing of 
electronic diffusion in disordered systems [I]- Being an 
interference wave phenomenon, this concept has been ex- 
tended to light, acoustic waves, and even Bose-Einstein- 
condensed matter waves Q. In the case of light, its 
vector character, together with the recent development 
of metamaterials with unusual electromagnetic proper- 
ties such as the possibility of optical magnetism, should 
lead to interesting particularities in AL. Indeed, the vec- 
tor character of light is at the origin of a polarization- 
induced anomalous delocalization effect in ID disordered 
systems lipl, for which the vast majority of states is lo- 
calized The presence of metamaterials in a ID su- 
perlattice also leads to light delocalization in the long- 



wavelength limit [12|. However, given that in Ref. [12 1 



the metamaterial was treated as non-dispersive (which 
leads to negative electromagnetic energy density isl. 14\) 
and that only normal incidence was considered, a more 
thorough investigation is clearly in order. 

The aim of the present Letter is to discuss light 
propagation through randomly perturbed ID photonic 
heterostructures composed of alternating layers of non- 
dispersive right-handed (RH) materials (labeled A) and 



dispersive left-handed (LH) metamaterials (label M) . We 
will establish that allowing for dispersion of both dielec- 
tric permittivity and magnetic permeability of the M lay- 
ers, together with oblique incidence leads to novel fea- 
tures of AL in ID. In particular, we find unexpected 
extended states at Brewster angles for both TM and 
TE configurations, a high transmission peak at the very 
edge of a bandgap, and light delocalization in the low- 
frequency regime at finite angles of incidence. 

We model our system as a stack of alternating layers 
of air (eA = y^A = 1) and of a Drude-like metamaterial, 
with responses for the dielectric permittivity and mag- 
netic permeability of the M layer given as [l5| 

eA/(t^) = eo ^, AiM(w) = Mo r, (1) 

such that Ve = uie/ {'^t^s/^o) and Vm — uJra/ {'^t^^/Jj^) are 
the frequencies associated with the electric and magnetic 
plasmon modes, respectively. We have followed previous 
work [15[ and use eo = 1.21 and //q = 1.0 in Eq. (1). 
Disorder is introduced by allowing the widths of the A 
and M components at the j-th layer to fluctuate around 
their respective mean values, a and b: aj = a -\- 8^ and 
bj ~ b + Sj^ , where the random variables S^'^' are homo- 
geneously distributed in the interval [—A/2, A/2]. The 
localization length f is calculated numerically using the 
standard definition jllj . 
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where T is the transmission coefficient, and L is the total 
stack length, L = + ^i); with N being the total 

• ) denotes configurational 



number of double layers; 
average. 
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The transmission coefficient is given by jl| 

2Z 



T 



(3) 
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where Z = cos 6, with 6 being the angle of incidence, and 
the elements of the transfer matrix Mij are defined as 
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with 



M 



i fx sin cosQxX 



, x = A, M, (5) 



Qx = {u}/c)uxiu},9),Uxiu},9) = ^ exiu})fix{(^) ~sin^6; for 
incident transverse electrical (TE) and transverse mag- 
netic (TM) waves, the coefficients fx are 
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For an infinitely periodic structure without disorder, Eqs. 
^ and (O lead to the dispersion relation 

cos{kd) = cos{qAa) cos(gj\/6) — sm{qAa) sin((7j\/6), 

(7) 

where F± = (fA/fM) ± (/a///a), with fx corresponding 
to TM or TE waves; d = a + 6 is the period of the system 
and k is the Bloch wavevector along the direction of the 
axis of the periodic photonic crystal. 

The numerical simulations are supplemented by a gen- 
eralization (to the case of oblique incidence) of an an- 
alytic expression for ^ derived by Izrailev and Makarov 
(IM) for bilayercd photonic structures [l^l , valid for weak 
disorder (small fluctuating widths in our case). 
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tac^Asin (qAffc) + gAfC^Af sin ilAa)], 

(8) 

where = (Sx)- For homogeneous random perturba- 
tions with the same amplitude on both layers, one has 
al=al = Ay 12.^ 

Let us first examine the band diagram for a perfectly 
periodic structure, as given by Eq. ([7]), and shown in Fig. 
[TJa) for TE waves. Dispersion leads to the appearance 
of two non-Bragg gaps: one at a frequency uj ~ An GHz, 
corresponding to the vanishing of the average refraction 
index of the structure, and the other, at w ~ Gtt GHz, 
which only occurs for oblique incidence, 0^0, and is 
a consequence of the excitation of electric and magnetic 
(for incident TM and TE waves, respectively) plasmon 
polaritons [l^l . Figure [TJa) also shows the usual Bragg 
gap at w ~ IStt GHz. Figure [ijb) shows the localiza- 
tion length ^ for TE waves as a function of frequency, 
for 6 = 7r/6. Before discussing the actual results, two 
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FIG. 1. (Color online) (a) Band diagram for TE waves, at 
incidence angle 6 = ti/Q, for a perfectly ordered superlattice, 
as given by Eq. (b) Localization length ^ (in units of 

the system size) for TE waves calculated numerically [Eq. 
0] for A'' = 5, 000 double layers and 100 realizations. The 
parameters used are: a = 6 = 12 mm, ujm = LUe ~ Gtt GHz, 
eo = 1.21, Ho = 1, and A = 0.5 mm (open red circles) or 
A = 12 mm (open blue squares). The red solid and blue 
dashed lines correspond to Eq. (O for A — 0.5 mm and 12 
mm, respectively. 



comments are in order: first, we have confirmed that the 
system is self-averaging, i.e., the behavior of ^ calculated 
from a single realization of disorder, for a system made 
up of a sufficiently large number of layers, does not differ 
significantly from that obtained through Eq. ([2]), consid- 
ering many disorder realizations for a system not so large; 
second, the overall agreement between the simulation re- 
sults and those from Eq. ^ is excellent for weak disorder 
(see data for A = 0.5 mm in Fig. [1]), but deviations oc- 
cur as disorder increases, as expected, as illustrated in 
Fig. [Ijb) for A = 12 mm. A comparison between Figs. 
[IJa) and (b) indicates that for weak disorder the dips in ^ 
correlate with Bragg and non-Bragg gaps alike; however, 
as the strength of disorder increases, the Bragg gap at 
a; ~ IStt GHz is washed out, leading to the disappear- 
ance of the associated dip in ^. The TM waves follow the 
same pattern. 

Most importantly. Fig. [TJb) shows that there are some 
specific frequencies at which ^ reaches anomalously high 
values, larger than the system size, hence leading to light 
dclocalization. This is a manifestation of the "Brewster 
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FIG. 2. (color online) The function g{uj), given by the RHS of 
Eq. (|9]), as a function of frequency for TE (full curve) and TM 
(dashed curve) waves. For a given angle of incidence, 9, the 
Brewster anomalies occur at frequencies satisfying cos^ = 
g{i->j)- As an example, the horizontal dotted line corresponds 
to cos^7r/6, whose intersections (indicated by arrows) with 
g{oj) yield the frequencies for Brewster anomalies for 9 = 7r/6. 
The media parameters are the same as in Fig. [1] 
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anomaly" [l^l, corresponding to the situation in which 
no reflection occurs in ID disordered optical systems at 
specific incident angles (Brewster angles, 9b)- To see 
how this comes about, one first considers non-magnetic 
materials {/im = 1); and a TM wave incident on an in- 
terface at an angle 9b'- there is no refiected component, 
since the induced electric field cannot radiate along its 
own axis. Multiple refiections are therefore eliminated 
along the medium, and localization is suppressed. For 
the metamatcrials considered here (/Um, cm 1)i the re- 
flected field is a result of radiation from both electric and 
magnetic dipoles, thus explaining the presence of sharp 
peaks in ^ both for TM and TE waves; Brewster anoma- 
lies for TM and TE waves have been observed at RH-LH 
interfaces [l^ . It is important to remark that the Brew- 
ster anomaly may occur at the vicinity of a band gap, 
so that the presence of disorder does not halt light prop- 
agation at band edges. Furthermore, at the band edge, 
group velocities may be very small and one could observe 
slow light propagation. 

Further insight into the frequency dependence of the 
Brewster angles can be obtained by examining the con- 
ditions for their occurrence within the framework of the 
generalized IM theory. According to Eq. (|5]), diverges 
when — > 0, which occurs for values of 9 such that 



cos 



9 - g{^) = J]—. — 



(9) 



where p{ui) = ij.m{^) for TE waves, and p(a;) = eA/('^) 
for TM waves. To illustrate this condition. Fig. [2] shows 
the frequency dependence of the Brewster angles, as cal- 
culated from Eq. ([5]), with the dispersive relations given 



FIG. 3. (color online) Localization length ^ (in units of the 
system size) for TE (a) and TM (b) waves as a function of the 
vacuum wavelength A, obtained from our numerical simula- 
tions, for different angles of incidence (solid lines, from top to 
bottom): — 0, tt/IOO, vr/12, n/6, and n/3. Vertical arrows 
locate the Brewster modes, and the dashed line corresponds 
to the asymptotic behavior ^ cx A® predicted in Ref. . The 
media parameters are the same as in Fig. [T] with A = 0.5 



by Eq. ([1]), for both TM and TE configurations. One sees 
that, as 9 increases, one goes from a regime in which the 
anomalies are present for three frequencies (correspond- 
ing to two TE and one TM mode, c.f. Fig.[2]for 9 = 7r/6) 
to one in which an additional TM mode appears. Equa- 
tion (jni) also puts in evidence that it is crucial to take 
dispersion into account in order to determine the values 
of Brewster angles in metamatcrials. 

Let us now discuss the behavior of the localization 
length over a wide range of frequencies, for different an- 
gles of incidence, 9. Figure 12] shows the simulation results 
for ^, for the specified set of parameters; from now on 
we refer to 27rc/a; as the (vacuum) wavelength A. First 
we note that, for normal incidence, ^ oc A^, similar to 
the behavior found in Ref. [l2j. However, this power- 
law behavior is already lost for very small angles of in- 
cidence (e.g., for 9 = tt/IOO); nevertheless, ^ may be 
larger than the system size in the long wavelength limit 
for small angles, 9 < 7r/12, thus suppressing localiza- 
tion. The long-wavelength delocalization for normal in- 
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cidcnce is due to the fact that we assumed = ojm- 
Indeed, the low-frequency Umit of the diverging condi- 
tion, Eq. (O, yields cos 6 ~ uj^l'^m for TE waves {ojml^^e 
for TM waves), so that = 0; by the same token, for 
We 7^ Wjn, light delocalization in the low- frequency regime 
occurs at 6* ~ cos^^[we/a;„i]'''^, which highlights the im- 
portance of allowing for oblique incidence. For larger 
angles {9 > 7r/12), the only mechanism for suppression 
of localization corresponds to the Brewster anomalies at 
well defined wavelengths: for TE modes and for a given 
angle of incidence, localization is suppressed at two dif- 
ferent wavelengths, one in the negative- and one in the 
positive-refraction regions of the metamaterial; for the 
TM modes, as discussed in connection with Fig. one 
must reach a critical angle of incidence before a Brewster 
anomaly emerges in the positive-refraction regime. 

Finally, as dispersive metamaterials are intrinsically 
absorptive, comments regarding losses are in order. The 
resonant nature of the electromagnetic response of the 
metamaterial, notably the magnetic one, is the main 
source of dissipation [2^; other sources of dissipation, 
such as surface roughness and interface effects, do ex- 
ist and constitute a major hurdle in the development of 
metamaterials. In the present Drude-like model, absorp- 
tion may be introduced phenomenologically through the 
replacement — >■ uj(uj -\-i^e,m)-: where 7e,m are the elec- 
tric and magnetic loss factors; we have found that the 
Brewster modes are indeed smeared out when absorp- 
tion effects are included. However, according to Eq. 
one may conceive realizing a metamaterial with dielectric 
and magnetic responses such that Brewster modes would 
show up in a spectral region in which losses are not too 
severe; in so doing, they could still lead to observation 
of light delocalization. Moreover, important progress in 
mitigating losses in metamaterials has been achieved re- 
cently, such as incorporation of gain 2l| , and the use of 
optical-parametric amplification pol 2^ . which has the 
advantage of being tunable in a wide negative index fre- 
quency range. These advances have led to significant 
attenuation and complete compensation of losses; even 
lasing in metamaterials has been observed There- 
fore, a combination of these loss-control mechanisms with 
the routes proposed here should lead to a wide range of 
possibilities for control of light fiow and Anderson local- 
ization in disordered metamaterials. 

In conclusion, wc have studied light propagation 
through ID disordered structures composed of alternat- 
ing layers of non-dispersive and dispersive metamateri- 
als; the disorder consists in allowing for randomness in 
the layer widths. We have found that the inclusion of 
dispersion, together with the possibility of oblique in- 
cidence, lead to novel routes for suppression of Ander- 
son localization of light, in addition to the long wave- 
length limit for normal incidence [l^ . First, for a given 
choice of dispersive metamaterials, oblique incidence at 



9 ~ cos~"'^[a;e/wm]^^ leads to delocalization in the long 
wavelength limit. The second route is provided by the 
Brewster anomalies, corresponding to a diverging local- 
ization length at well defined frequencies. Interestingly, 
these anomalies can occur for frequencies in ranges such 
that the metamaterial exhibits either negative- or posi- 
tive refraction; thus, negative refraction is not a neces- 
sary condition for delocalization, but the use of dispersive 
materials is. These results indicate that control and fine- 
tuning of multiple scattering of light could be achieved 
by a careful combination of angle of incidence and a ju- 
dicious choice of metamaterials. 
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